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We investigate ferromagnetic ordering in an itinerant ultracold atomic Fermi gas with repulsive 
interactions and population imbalance. In a spatially uniform system, we show that at zero tem- 
perature the transition to the itinerant magnetic phase transforms from first to second order with 
increasing population imbalance. Drawing on these results, we elucidate the phases present in a 
trapped geometry, finding three characteristic types of behavior with changing population imbal- 
ance. Finally, we outline the potential experimental implications of the findings. 



PACS numbers: 03.75.Ss, 71.10.Ca, 67.85.-d 



I. INTRODUCTION 

Feshbach resonance phenomena provide unprecedented 
control of pair interactions in degenerate atomic Fermi 
gases P, [|| . This feature has allowed extensive studies of 
pairing phenomena in two-component Fermi gases pro- 
viding access to the crossover between a Bose-Einstein 
condensate (BEC) of molecules and the Bardeen-Cooper- 
Schrieffer (BCS) state of Cooper pairs @, 0, S, @|. Al- 
though the emphasis of experimental investigations has 
been primarily on the problem of resonance superfluidity, 
interacting Fermi gases support other strongly-correlated 
phases including itinerant ferromagnetism. 

In solid state condensed matter systems, the problem 
of itinerant ferromagnetism has a long history dating 
back to the pioneering studies by Stoner 0] and Wohl- 
farth [8J. These early investigations proposed that, at 
low enough temperatures, a Fermi gas subject to a re- 
pulsive interaction potential could undergo a continuous 
phase transition into an itinerant spin polarized phase Q . 
This Stoner transition reflects the shifting balance be- 
tween the potential energy gained in spin polarization 
through Pauli exclusion statistics, and the associated cost 
in kinetic energy. Subsequent studies showed that fluctu- 
ations in the magnetization at low temperatures drive the 
second order transition first order at low enough temper- 
atures [ISIil EUS Ellil- Such behavior is born out 
around quantum criticality in a variety of experimental 
solid state systems including ZrZn? [la . [T3] , UGe? fisj ] , 
MnSi [Ilj2p,[2l|,[23,[ll], CoS 2 [if, YbRh 2 Si 2 df, and 
SrRu03 [26|] . When subject to a magnetic field, the at- 
tendant increase in Zeeman energy results in the bifurca- 
tion of the tricritical point separating the region of first 
and second order ferromagnetic transitions into two lines 
of mctamagnctic critical points. 

In the following, we will explore the potential im- 
plications of this itinerant magnetic phase behavior on 
the equilibrium properties of strongly interacting two- 
component atomic Fermi gases; here we refer to the 
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pseudo-spin associated with the hyperfme states char- 
acterizing the two atomic populations. However, in con- 
trast to the solid state system, the application of these 
ideas to the atomic Fermi gas must address the features 
imposed by the trap geometry, and the constraints re- 
sulting from the inability of particles to transfer between 
different spin states [27]]. As a result, in the general case, 
one must consider atomic Fermi mixtures in which an 
effective spin polarization is imposed by population im- 
balance [28l . |29(. The potential for itinerant ferromag- 
netism in atomic Fermi gases has been already addressed 
in the literature, Zhang et al. [HJ showed that an opti- 
cal lattice would reduce the repulsive interaction strength 
required to realize ferromagnetism, and Sogo and Yabu 
[30( studied a trapped system in the Thomas-Fermi ap- 
proximation. Subsequently, Duine and MacDonald [3 11 ] 
developed a diagrammatic pcrturbative expansion in in- 
teraction strength to address the phase behavior of the 
balanced two-component Fermi system. In the following, 
we will develop a functional integral formulation to ex- 
plore the phase behavior of the general population imbal- 
anced system. As well as providing access to the mean- 
field phase behavior of the system, such an approach al- 
lows for future considerations of the collective low energy 
spin dynamics of the spin polarized phase. Moreover, the 
theory provides a platform to explore the potential for 
the development of an equilibrium spin textured phase 
recently conjectured in relation to the solid state sys- 
tem [13, m, s m, hi . 

The paper is organized as follows: In Sec. |TT] we de- 
rive an expression for the thermodynamic potential of 
the system as a function of the local density and in-plane 
magnetization fields. To address the important effects of 
spin-wave fluctuations on the nature of the equilibrium 
phase diagram, we will explore the rcnormalization of the 
mean-field equations keeping those terms that are second 
order in the coupling strength, g. Using this result, in 
Sec. IIII Al we analyze the phase diagram of the spatially 
uniform system as a function of the interaction strength, 
g, and chemical potential shift. Finally, in Sec. IIII 51 we 
explore in detail the phase behavior of the magnetic sys- 
tem in the atomic trap geometry. 
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II. FIELD INTEGRAL FORMULATION 

Expressed as a coherent state path integral, the quan- 
tum partition function of a population imbalanced two- 
component Fermi gas is given by 

Z = Jvipcxp — J drdr ^ Vv(— id T + | — <TAp)ip a 



<7={T,!} 



drdr gip^tp tip ^ 



(1) 



where ip a {T, r) and ip a {T, r ) denote Grassmann fields, (3 = 
l/fceT is the inverse temperature, and £ = p 2 /2m — p. 
Here we have used a pseudo-spin index, a G {1,4}, to 
discriminate the two components. As independent par- 
ticles (with no interconversion) , the density of the two 
majority/minority degrees of freedom must be specified 
by two chemical potentials. For convenience, it is helpful 
to separate the chemical potentials into their sum and 
difference; p + Ap for up-spin and p — Ap for down-spin. 
In this representation, population imbalance may be ad- 
justed through the chemical potential shift, Ap. Note 
that, although population imbalance is synonymous with 
a global pseudo-spin magnetization, a spontaneous sym- 
metry breaking into an itinerant ferromagnetic phase can 
still develop with the appearance of a non-zero in-plane 
component of the magnetization. Finally, we suppose 
that the strength of the repulsive s-wave contact interac- 
tion, <?<5 3 (r), can be tuned using a Feshbach resonance. 



A. Hubbard-Stratonovich decoupling 

To develop an effective low-energy theory for the Fermi 
gas, it is convenient to decouple the quartic contact inter- 
action by introducing auxiliary bosonic fields, p and cp, 
conjugate to the local density 5Z a ={| j.} V'qV'q an d mag- 
netization J2 a /3={T 1} i^uVapipp respectively, setting 



Z = / V^VpVipexpi - / drdr 



l{<t> 2 



r) 



(G 1 + gp)S a( 3 - (Ape z + gcj>) ■ cr a p 



(2) 



Here Go = (— i5* + £) _1 defines the Green's function 
of the non-interacting system, and er denotes the vec- 
tor of Pauli spin matrices. Note that, without de- 
coupling in both the Hartree and Fock channels, one 
would subsequently encounter unphysical diagrammatic 
contributions to the perturbative scheme developed be- 
low [37], [H, [39|. It is also the simplest approach that 
maintains spin rotational invariance of the Hamilto- 
nian, and leads to the correct set of Hartrcc-Fock equa- 
tions [i^, 5H- Then, integrating over the Fermi fields, 



one obtains the expression 

/ P0Pp2?Ve" /drdr9(02 ~ p2) 



Tr In (Gq 1 + gp - er • (Ape z + gep)) . (3) 



x exp 



At this stage the analysis is exact, but to proceed further 
one must employ an approximation. To orient our discus- 
sion and make contact with conventional Stoner theory, 
let us first consider a direct saddle-point approximation 
scheme. 



B. Stoner mean-field theory 

As well as the "effective" magnetization imposed by 
population imbalance, we anticipate the development of 
a spontaneous magnetization which will drive the axis of 
quantization away from the z-axis. We re-orient the axis 
of quantization to lie parallel to the net magnetization, 
denoted in mean-field theory (with over-bars) cp = cp ± + 
4> z e z , <p^ = (0 x ,(/) y ), and with this definition, the total 
magnetization of the system is given by M = Ape z / g+<p. 
Separately varying the action with respect to tp±_ an d 4> z 
one obtains, respectively, the saddle-point equations, 



(pV)- 1 Tr{G+-G-) 



(.90j 



(<70 z + Am) 5 



Ap 



where G ± = G 



g~p =F |A^e z + g<f>\, and V denotes 
the total volume of the system. Together, these equations 
admit two possible solutions: 

[cp ± = and M = <f> z ] : The total magnetization of the 
system can be ascribed to population imbal- 
ance with no spontaneous magnetization in-plane. 
Within this solution, M is a function of \g<t> z + Ap\, 
so it can be used to infer the chemical potential 
shift, Ap. 

[cp ± 7^ 0]: The total magnetization takes the form M = 

(cp ± + t^) 1 / 2 . Along z-axis, the magnetization is 
fixed due to population imbalance, with the ad- 
ditional magnetization developing within the x-y 
plane. In this case, the saddle-point solution trans- 
lates to the condition Ap = 0, i.e. no chemical 
potential shift is required to recover the fixed z- 
component of the magnetization due to the popu- 
lation imbalance; it is simply given by the resolved 
component of the total magnetization. 

The total population N = iV-f + N± can in turn be ob- 
tained from the variation 5S/S~p = 0. 

Expanding the action in interaction strength, g, S = 
g</>zTr(l + gG a G a ) = g$ 2 z {\ — gv), and one can ex- 
tract the familiar Stoner criterion [Z^, 51] for a pop- 
ulation balanced system, with v being the density of 



3 



states. For gv < 1 the state is unmagnetized, M = 0, 
and chemical potentials of the two Fermi surfaces re- 
main eq ual. If gv > 1 t hen the state is magnetized with 
M = yj {gv — l)/g 3 v". We also note that the Stoner 
criterion can be reformulated to account for population 
imbalance giving S = g(j> 2 (l — gv) — j A/i , leading to 
a transition at the same value of interaction strength as 
for the balanced system. Although, at this order, the 
saddle-point approximation predicts a continuous transi- 
tion to a ferromagnetic phase for the balanced system, it 
is well-established that fluctuations of the magnetization 
field drive the transition first order at low temperature 
[44] ]. This effect can be captured by retaining fluctua- 
tion contributions to second order in the interaction. In 
the following, we will explore the impact of fluctuations 
on the equations of motion associated with the uniform 
mean-field. 



C. Integrating out auxiliary field fluctuations 

To implement this program, it is convenient to param- 
eterize the Hubbard-Stratonovich fields into some, as yet 
undetermined, stationary (spatially uniform) values <fio 
and po, and fluctuations around them, cpa and pa- In- 
tegrating out these fluctuations, the goal is to obtain 
the renormalized mean-field equations for <fio and po re- 
taining contributions to second order in g. Substituting 
<f> = <fio + 0fl and p = po + pn into Eq. ©, and ro- 
tating the z-axis from the quantization direction to lie 
along the direction of uniform magnetization using the 
constant matrix T, one obtains 



Z 



VpfiVcpa exp 



drdr g((j>l - pi) 



x exp[TrlnG- 1 +Trln(]-|-stGT- 1 (l / 9 fl + o--0 fl )T)] , 

where now Gjjj 1 = Gq 1 + gpo =F \Ape z + g4>o\ denotes 
the elements of the inverse Green's function of the sys- 
tem at the level of the renormalized mean-field, G = 
diag(G+, G_). Then, expanding the action to second or- 
der in fluctuations, /Ofl(r, r) and 0fl(r, r), and performing 
the functional integral, one obtains the thermodynamic 
grand potential from the quantum partition function us- 
ing <I>g = ~/3~ l InZ, 



t 



TrlnG: 



TrlnGI^gi 



Po) 



1 



Trln(l-.g 2 II ++ n_ 



1 



+ - Trln (1 + gIl + _ + gU_ + + .g 2 n + _II_ + ) , (4) 

a result that is independent of the transformation T. 
Here we have defined the spin-dependent polarization op- 



erator, 



IW (w, q) = V G s (w', k)G s , (uj' - u, k - q) , 

UJ ,K 



where the sum on lu' runs over fermionic Matsubara fre- 
quencies. The term labeled (f) simply represents the 
thermodynamic potential of a non-interacting Fermi gas 
with shifted chemical potentials. The term labeled (_L) is 
due to transverse fluctuations of the magnetization field 
and coincides with that obtained in Ref. [45| . By con- 
trast, the term labeled (|j), corresponding to longitudinal 
fluctuations, differs from that obtained in Ref. [45| by the 
additional contributions from density fluctuation effects. 

To proceed, we now expand the potential <I>g to second 
order in g and perform the summations over Matsubara 
frequencies. Rearranging the momenta summations, one 
obtains 



*G = —5y J2 ln(l + C-^ ek ^) 



k 



+ g +gN+N- 




2g 2 x ^ n + (e k l )n_(e k2 )(l-n + (e k3 ))(l-n_(e k4 )) 
V ^ 



£ki + £k 2 — e k 3 — e k 4 



(5) 



where p s = p - gp + s|A^e z + gcp \, n s (e) = (1 + 
cxp(— f3(e — /z«))) , and N s = X)k n s(ek)- Conservation 
of momentum requires that ki + k 2 = ks + k4. Physi- 
cally, the numerator of the second order term indicates 
that the matrix element associated with the transition 
(ki, k 2 ) — > (k3, k.4) is proportional to the probability that 
states ki and k2 are occupied, whilst states k3 and k4 
are unoccupied. Following Pathria [46| (and the earlier 
discussion of Abrikosov and Khalatnikov [3), to renor- 
malize the unphysical divergence of the term in n 2 (e), 
labeled (0) close to resonance, we regularize the effective 
interaction at second order in scattering length a, 



5(ki,k 2 ) 



2fcpa 



n 2 v 2 V 2 ^ e kl + e k2 

k3,4 



e k 3 



e k 4 



where v = y/Jt/ v2tt 2 and kp = \/2mp. In a popu- 
lation imbalanced system the definition for the chemi- 
cal potential is that which gives the same total number 
of particles in the population balanced system, that is 
fcp = ^/37r 2 (n| + raj), where w\ and are the number 
of up and down-spin particles; this definition holds true in 
both the canonical and grand canonical ensembles. This 
regularization of the contact interaction exactly cancels 
the divergent terms in n 2 (e), labeled (<))• Furthermore, 
the terms in n 4 (e) are zero by symmetry. Finally, making 
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use of the symmetry in k3 and k4, one obtains From the thermodynamic potential we can compute the 

free energy per unit volume F = $g + So-={t !}(/■* + 

$ G = In ( 1 + c _ ^( ek_AIs ' j <jAn)N a . To consolidate terms entering the free energy 

P k we switch from the population imbalance pseudo-spin ba- 

J5 ={+'~} sis to the magnetization basis, retain contributions to or- 

Hk^a 2 ^ n +( e ki)^-(ek 2 ) (n+(e k3 ) + ?i_(e k3 )) 



2fc F q , 2 _ 2\ + 2fc F q ^ ^ dcr C((fc F a) 2 ), recall that if Ap = then M ^ 0, whereas 

^ ' txv + if A/z ^ then M = 0, and affect the rearrangement 



Kl,2,3 



(6) 



2fc F a 



A/ie z 



2fc P a 



AL 



- Pi) + ( 2J *rPo + Ap - |A/xe z 
A/i+lA^ + ^^ol)^- 



2kpa 



-Po 



Ape z 



2kya 



: 4>o 



AL 



Then, if we set 4>q = 4> + A<p and po = P + Ap, an ex- 
pansion in A<fi and Ap shows that the terms labeled (0) 
sum to zero to the accuracy of the free energy, 0((k F a) 2 ). 
Retaining the remaining contribution, the free energy re- 
duces to the form, 



F-- 



L 5>(l+e-«*-*>) + 5>* 



pv 



k 

{ + ■-} 



2kya 



-N+N- 



*={+,-} 



8fc|a 2 n + (e kl )n_(e k2 ) (n + (e k3 ) + n_(e k3 )) 



^2^3 ^ 

kl,2,3 



£ki + e k2 — e k 3 — e k4 



This expression coincides [55j with that obtained in 
Ref. [3l|. The method employed in the numerical calcu- 
lation of the summation over three momenta is described 
in App. El 



D. Magnetization 



To minimize the free energy and obtain the net mag- 
netization it is convenient to take the expression for the 
thermodynamic potential ^ and affect the shift of the 
field $ z i— > <f> z — Afiiris /2kp a . As a result, the thermody- 



namic potential takes the form 



= - — y lnfi + e-^*-"') 



2k F a 



4>a 



A/ie z 7w 



2h 



2kpa P 
Po 



2fcpa 



tt;/ 



n+(e kl )ri-(e k2 ) (n+(e k3 ) + n_(e k3 )) 

kl 2.3 



e ki + e k2 



e k 3 



e k 4 



(7) 



where, in response to the shift of 4> z , the factors of fj, s = 
H — 2kpapo/irv + 2kpas\4> \/ni' entering the definitions 
of N± and n± are now independent of Ap. The thermo- 
dynamic potential can be rewritten in terms of a func- 
tion of just the auxiliary fields and the chemical potential 
shift as $g = F(\4>o\) + 2kpa\(f)o - Ape z np /2fcF Q | 2 / 'ttv - 
2k~pap\j'Kv. 

In the grand canonical ensemble, the thermodynamic 
potential must be minimized with respect to the compo- 
nents of the auxiliary field giving 



4fcpa 



4>j 



= o, 



|0o| V ^ z ) V ^ z _ A/i7ri//2fc Fa 

(8) 

where <j>± = {<f>x,<f>y) so (po = 0j_ + 4> z e z . Following 
Sec. Ill Bl one may now identify the magnetization with 
the field <po — 7wApe z /2fcpa. If Ap = 0, then the system 
of equations is solved by either F' (\4>q\) / \ 4>o\+4k-pa/ nis = 
(the direction of spontaneous ferromagnetism in-plane 
remains undetermined), or 0q = 0. If A/i =/= then 
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FIG. 1: (a) shows the magnetization M as a function of 
population imbalance, P and interaction strength kpa — 
^/37r 2 (ri| + ni) m the canonical ensemble at T — at fixed 
species populations. The thick line traces system variation 
at P = which corresponds to trap profile {P/N = 0) in 
Fig- El (b) shows the phase boundary between "unmagne- 
tized" (UnM) and partially magnetized (PM) region and the 
line of saturation before the fully-magnetized (FM) region. 
Note that, by unmagnetized, we refer to the not in-plane mag- 
netization. The solid line denotes first order transitions, the 
dashed second order and saturation. 



4>±_ = 0, and the magnetization is set by the equation 
F'((f> z ) = 2(A/i — 2kpa4> z /TTi/) and is oriented along the 
z-axis. This behavior is analogous to what we saw in the 
mean-field analysis in Sec. Ill Bl Finally, as a consistency 
check, one may note that the expected degree of popula- 
tion imbalance can be recovered from the grand potential 
M = -d$ G /dAfi\ TV>N . 



is no net magnetization. As the interaction strength is 
increased, at kpa w 1.05 there is a first order phase tran- 
sition into a magnetized phase with M/N w 0.6. As 
fcpa is increased further the magnetization rises until it 
is saturated at kpa ~ 1.11. 

With increasing population imbalance, P, at kpa < 
1.05, where it is not energetically favorable for a spon- 
taneous magnetization to develop, the magnetization is 
forced to stay pinned to the minimum value set by 
the imbalance. With increasing interaction strength, at 
kpa « 1.05 there is a first order transition and the magne- 
tization jumps to M/N ss 0.6. This feature is consistent 
with the findings of the Stoner mean-field theory that 
the transition interaction strength found is independent 
of population imbalance. If the population imbalance is 
greater than P/N > 0.6 then the magnetization takes the 
value of the spontaneous magnetization projected onto 
the sheet of minimum magnetization caused by the pop- 
ulation imbalance. 

From these results, one can infer the corresponding 
zero temperature phase diagram Fig. [TJb) . Characteriz- 
ing the phase behavior by the strength of the in-plane 
magnetization and the degree of polarization, the phase 
diagram divides into three distinct regions. At low inter- 
action strength the system is not spontaneously unmag- 
netized, though there can be a magnetization fixed by 
the population imbalance. Then, at increased interaction 
strength the system become partially magnetized either 
through a first order (at low population imbalance) or 
a second order phase transition. At interaction strength 
above kpa > 1.11 the magnetization saturates. 

To address the properties of the population imbalanced 
system in the grand canonical regime, we will divide our 
discussion between the uniform and trap geometries. In 
Scc. lIII Al we will address the properties of a uniform sys- 
tem where the chemical potential /i and shift Afi are held 
constant (allowing the species populations to effectively 
interchange). Drawing on these results, we will then dis- 
cuss the phase behavior in a harmonic trap in Sec. IIII Bl 



III. POPULATION IMBALANCE 



A. Uniform system 



With the formal part of the analysis complete, we will 
now apply these results to explore the implications of fer- 
romagnetism in the atomic Fermi gas. To begin, let us 
consider the phase behavior of the system in the canoni- 
cal ensemble working at fixed particle number. The vari- 
ation of the total magnetization, \M\, as a function of 
interaction strength and particle imbalance can be found 
by minimizing the free energy at fixed particle number. 
The results are shown in Fig. [TJ To ensure that the free 
energy is locally minimized rather than just being at a 
stationary value [47] , the curvature was examined nu- 
merically. In the balanced Fermi gas, P = 0, the results 
shown in Fig. (Ha) recapitulate those discussed by Duine 
and MacDonald [3l|. In particular at zero temperature, 
when the interaction strength is small, kpa < 1.05, there 



In the spatially uniform system, when the chemical po- 
tentials of the two species are fixed, for each value of the 
interaction strength kpa and relative shift in chemical 
potential Afi/n, from the free energy one can obtain the 
phase corresponding to minimal thermodynamic poten- 
tial. Applying this procedure, the resulting phase behav- 
ior is shown in Fig. O For A/i//i = and small interac- 
tion strength kpa < 1.05 there is no magnetization. As 
the interaction strength is increased, at kpa ks 1.05 in 
the canonical regime Fig. [1] there is a first order phase 
transition into a fully-magnetized state. Working at 
fixed chemical potential [Fig. [U^a)], the phase transition 
straight into a saturated state increases the number of 
particles, which in turn increases the effective interaction 
strength to kpa w 1.25 (calculated using the chemical po- 
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(P/N = 0.8) 



{P/N = 0.4) 
■(P/N = 0) 




(P/N = 0.8) 



(P/N = 0.4) 
(P/N = 0) 



FIG. 2: (Color online) (a) shows the variation of magneti- 
zation M as a function of chemical potential shift Ajj, and 
interaction strength k^a in the grand canonical ensemble at 
T = 0. The thick lines correspond to trap profiles at P/N = 0, 
P/N = 0.4, and P/N = 0.8 in Fig. [3] In the region where 
magnetization is undefined there is phase separation. The 
lower set of diagrams show the phase boundaries (and satu- 
ration line) between "unmagnetized" (UnM), partially mag- 
netized (PM) and fully-magnetized (FM) regions, as well as 
the region of phase separation. 



tcntial for a non-interacting system with the same total 
number of particles). This leads to an intermediate re- 
gion of phase separation in the grand canonical regime. 
At kpa < 1.05 as the chemical potential shift is increased 
up to A/i/yU = 1, the magnetization increases up to its 
maximum saturated value as the Fermi surfaces become 
more unbalanced. At A/i//i > 1 the chemical potential of 
the minority spin species is negative so only the majority 
spin species remain and the system is fully magnetized. 
With a chemical potential shift the region of phase sep- 
aration corresponds to the first order phase transition in 
Fig. [T] The corresponding phase diagram showing the 
regime of two- phase coexistence is shown in Fig. EJb). 

Finally, if the system has an imposed density and pop- 
ulation imbalance, and the chemical potentials are free to 
vary, then there are two possibilities: Firstly, the sponta- 
neous ferromagnetism is sufficient to provide the popula- 
tion imbalance and any excess magnetization lies in the 
plane. This corresponds to a point on the line A/i = 
in Fig. [2J The second possibility is that spontaneous 
ferromagnetism is not sufficient, and so there is an addi- 
tional chemical potential shift A/t ^ 0. In this case the 
magnetization then points along the direction of popu- 
lation imbalance. This is consistent with the findings in 
Sec. Ill Dl For a given interaction strength, the magne- 



tization increases with chemical potential shift to satu- 
ration, so there is always a chemical potential shift that 
will give a suitable population imbalance. 



B. Trapped system 

Using the insight gained from the study of the uni- 
form system, we can now explore an atomic Fermi gas 
in the physical system — a potential trap. Without 
loss of generality we take f (J.) to represent the ma- 
jority (minority) species of atoms. We focus on a har- 
monic trap, with rescaled spatial coordinates to ensure 
a spherically symmetric trapping potential, V(r) ~ r 2 . 
Furthermore, we make use of the local density approx- 
imation in which the chemical potential of both species 
A t eff,a'(r) = ^o.ct — ^( r ) are renormalized by the same 
trapping potential. Although there is some experimen- 
tal evidence [H, [49 |_tliat the local density approximation 
might not be valid 10,[5lJ in some setups, we believe that 
its application here will correctly address the qualitative 
phase structure. The chemical potentials are regarded to 
be locally fixed, therefore the local phase is that of the 
uniform system in the grand canonical regime examined 
in Sec. lIII Al With a constant chemical potential shift A/i 
and interaction strength <?, but varying effective chemical 
potential /x, the system follows the trajectory kpa oc JJi 
and Afi/fi oc in the grand canonical regime shown 
in Fig. [21 If the chemical potential is large, the system 
spontaneously becomes ferromagnetic, and the magneti- 
zation is saturated; if the chemical potential is small, the 
relative chemical potential shift is large ensuring the mag- 
netization is again near saturation. The locus in Fig. [2] 
shows that, in the intermediate region, the magnetiza- 
tion can develop a minimum depending on the degree of 
population imbalance. 

To understand the behavior in the trap geometry, one 
should note the following: If the degree of equilibrium 
pseudo-spin magnetization is in excess of that imposed by 
total population imbalance alone, the analysis of Sec. Ill Dl 
tells us that some component of the spontaneous magne- 
tization lies along the z-axis with the remainder oriented 
in the x-y plane. If, however, net population imbalance is 
large, then A/i ^ and no in-plane magnetization devel- 
ops. Here one may identify three characteristic behaviors 
with radial density profiles shown in Fig. [3J The first 
{P/N = 0) has in-plane magnetization, and the others 
do not. The second {P/N = 0.4) has a first order tran- 
sition and non-zero phase separation whereas the third 
{P/N = 0.8) is always fully magnetized due to strong 
interactions. The three plots all have the same central 
chemical potential. 

The first possibility shown in Fig. ^P/N = 0) is at 
small population imbalance, involving the development 
of a spontaneous magnetization which is in excess of what 
can be absorbed by population imbalance alone, in this 
case A/i = and some magnetization lies in the plane. 
At small radii, where the interaction strength kpa > 1.25 
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FIG. 3: The density of particles at radius r in a trap potential 
profile at three different values of total population imbalance. 
The variation in the local particle density N is shown by the 
solid line, and the variation in the local magnetization M is 
shown by the dashed line. The plot densities are renormal- 
ized by their central density, No, radii by the outer radius, 
ro, of non- interacting particles with the same average inner 
chemical potential, fj,o. The upper panel shows the effective 
k^a in the P/N = case. 



is greater than the limit for ferromagnetism, the results 
of the uniform system (Sec. IIII Ap show that there is sat- 
urated ferromagnetism in the plane and a normal compo- 
nent that provides the fixed population imbalance. Fol- 
lowing this there is a region of phase separation and then 
at k-pa w 1.05 there are equal particle densities and no 
magnetization. The outer edge of the particle distribu- 
tion of both species is where /i = V(r Q ). 

In the second scenario shown in Fig. ^P/N = 0.4) 
the spontaneous magnetization is not sufficient to pro- 
vide population imbalance alone, in this case we require 
Afi =/= 0, and all magnetization is oriented along the axis 



of population imbalance. From the trap center the pop- 
ulation imbalance is first fully saturated, followed by a 
region of phase separation, into a region of partial magne- 
tization. This causes the minority spin particles to have 
a sharp maximum number density at r/ro as 0.6, and the 
magnetization to have a corresponding minimum; this 
counters the intuitive expectation that number density 
should rise towards the trap center due to the increas- 
ing effective chemical potential. As the effective chemi- 
cal potential continues to fall with increasing radius, the 
minority spin species population falls more rapidly than 
the majority and magnetization increases. At a large 
radius, the chemical potential of the minority spin par- 
ticles reaches zero before the majority spin so there is a 
thin shell containing only majority spin particles at the 
outside and so is fully magnetized. 

The third possibility shown in Fig. gfP/iV = 0.8) is 
that the locus in Fig.[2]does not cross the first order tran- 
sition and region of phase separation. At A.fj,/fx < 1 the 
system is fully magnetized due to the strong interactions 
between particles. At Ayu/yU > 1 the system is fully mag- 
netized due to there being no minority spin particles. In 
the intermediate regime there is a narrow band where the 
system is partially polarized. The majority spin species 
exists out to greater radius than in cases (P/N = 0) and 
(P/N = 0.4) because A/j is larger so a greater potential 
at a larger radius is required to give the majority spin 
species zero effective chemical potential. 



IV. DISCUSSION 

To conclude, let us now consider four methods of 
how spin magnetization could be detected experimen- 
tally. Firstly, the interaction energy can be estimated 
by studying the expansion of the gas [52] . Time of flight 
measurements of the expanding cloud with no external 
magnetic field B = are ballistic and so can provide 
the initial kinetic energy. If the magnetic field is present, 
B =/= 0, then interactions are significant during the expan- 
sion. Collisions ensure that all of the interaction energy is 
converted into kinetic energy so the measurements reflect 
the total released energy. Taking the difference between 
the B ^ and B = measurements therefore probes the 
interaction energy. An unmagnetized gas has interaction 
energy whereas the fully magnetized gas has zero interac- 
tion energy so time of flight measurements should allow 
the ferromagnetic state to be detected. 

Radio frequency spectroscopy [53| allows one to probe 
the spatial variations of scattering lengths by exciting 
the atoms from one spin state |1) into some other state 
1 3} whilst leaving the atoms in the second spin state |2) 
unaffected. The presence of atoms in state |2) shifts the 
resonance U13 by A^i3 = 2712(023 — 021), where is 
the scattering length between states \i) and \j). Mea- 
surement of the resonance shift could allow the spatial 
distribution of the individual species to be probed. The 
presence of the ferromagnetic state could be inferred by 



8 



looking for the characteristic density profiles outlined in 

Sec. nns 

A third simple method of detecting a ferromagnetic 
transition could be to monitor the size of the atomic 
cloud. In a harmonic trap the cloud size is proportional 
to the square root of the Fermi energy. Therefore, the 
size of the fully-magnetized state is 2 1 / 3 larger than the 
unmagnetized. 

On the repulsive side of the Feshbach resonance three- 
body collisions can result in the formation of a molecular 
bound state of two atoms that might destroy the atomic 
gas before it has time to undergo ferromagnetic ordering. 
To overcome this obstacle an atomic gas spin could be po- 
larized along the magnetic field direction and an RF w/ 2 
pulse applied to rotate all the spins into the plane [531 ] - 
The rate of precession of the spins is set by the magnetic 
field strength, which varies across the atomic gas due to 
field inhomogeneitics. The precession rate of the atoms 
would however be kept locked together by the ferromag- 
netic interaction. Furthermore the ferromagnetic phase 
has an antisymmetric wave function which inhibits colli- 
sions and so prevents the formation of molecular bound 
states. A signature of ferromagnetism is therefore the 
absence of molecular bound state formation. 

We now outline two possible ways to further our anal- 
ysis. The first order phase transition leads to discontinu- 
ities in the density and magnetization leading to phase 
separation. Such behavior could lead to a breakdown of 
the local density approximation, a potential source of in- 
accuracy in our analysis. This could be fixed through 
inclusion of a surface energy. 

The second is to investigate the possibility that mag- 
netic texture could develop. Textured modes may have 
been seen via the possible formation of a CDW/SDW 
in experimental results on the analogous solid state sys- 
tems of itinerant electron ferromagnets UGe2 [HI, |34| . 
Ca 3 Ru207 [54[, and MnSi Our general formalism 

should be able to be extended to include the possibility 
of a textured phase which lies beyond the first order line 
in the putative paramagnetic regime. 

In conclusion we have developed a general formalism 
to describe itinerant ferromagnetic transitions in two- 
component fermionic cold atom systems with repulsive 
interactions, and potential population imbalance. At 
low population imbalance, we predict that the first order 
transition that characterizes the balanced system per- 
sists. However, when the imbalance is large the transi- 
tion becomes continuous. In the trap geometry we found 
the first order phase transition led to discontinuities in 
density and magnetization. Up to a critical total popu- 
lation imbalance, set by the possible total magnetization 
following a first order transition, the phases in the trap 
had the same density and magnetization profiles with in- 
creasing population imbalance, but in-plane magnetiza- 
tion fell. With population imbalance above this level, the 
system requires a chemical potential shift to generate a 
population imbalance; however there is still a small range 
over which a first order phase transition is seen. In the 




FIG. 4: The re-parameterization of momenta used to en- 
sure momentum conservation. ki,2,3,4 represent the momenta 
appearing in the original integral, whose separate sums are 
qi2 = ki + k2 and q34 = kj + k4, which the Dirac delta func- 
tion in Eq. (|A1|I will ensure that qi2 = q34. 612 represents 
the angle between ki and k2, k^ is the vector perpendicular 
from qi2 to ki and k2, and k^ is similarly defined. 

two latter cases the local population imbalance displayed 
a characteristic minimum with radius. 
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APPENDIX A: COMPUTATIONAL ANALYSIS 
OF MOMENTUM SPACE INTEGRAL 

An important integral Eq. (J7]) encountered in this pa- 
per has the form 

F(|ki|, |k 2 |, |k 3 |, M)d(ki + k 2 - k 3 - k 4 ) 
xdkidk 2 dk 3 dk4 . (Al) 

To evaluate this integral one could substitute k 4 = 
ki + k 2 — k 3 , and then integrate over three parameters 
representing the lengths of vectors ki, k 2 , and k 3 , and a 
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minimum of three relative angles between these vectors, 
giving a total of six integration parameters. However, 
numerical integration generally becomes more prohibitive 
with increasing number of dimensions. Since the function 
F depends only on the magnitude of the momentum, the 
scheme outlined below allows us to perform the angular 
integration separately of the function and leave a numer- 
ical integral over just the four dimensions of the vector 
lengths. 

The integral is re-parameterized according to Fig. 31 
qi2 = k!+k2 and = k 3 +k4, the vector perpendicular 
from qi2 to ki and k2 is k^, and k^ 4 is similarly defined. 
The vector k^ 2 has length given by 

kh = Y q \J^ 2 {kl + kl)- q ^-{kl-klf . ( A2) 

We first concentrate on calculating the angular compo- 
nent just of the integral over ki and k2, the angle be- 
tween these vectors is #12. The phase space volume of 
the angular integral is 



which is the expected result. A similar procedure is used 
to parameterize the separate integral over the angular 
components of k 3 and k 4 into q 3 4. 



The original integral Eq. (|Alj) is now re-written in 
terms of the parameters qi2 and 534 using Eq. (jA4|. Mo- 
mentum conservation is required by the presence of the 
Dirac delta function (5(ki + k 2 — k 3 — k 4 ) = 5(qi2 — q34), 
however the q\2 and q 3 4 parameters introduced are just 
scalar quantities. The momentum conservation require- 
ment is implemented by demanding that the two scalar 
integration parameters are equal, which sets the two in- 
tegration parameters equal, q\2 = 9 3 4 = q so there is just 
one integral over parameter q remaining. However, this 
introduces an extra angular degree of freedom (the angle 
between qi2 and q 3 4). In order to compensate the inte- 
grand is divided by the extra phase space volume of the 
angular integration, Anq 2 . Wc then obtain 



sin0i2<i0i2 



£-L / 1 1? 2 — h 1 - 2 

ft 12 / r, 1 / 1 12 




h.2 _ h.±2 \ 
2 12 W 



h.2 _ L,±2 1 \ h.2 _ jL_L2 
ft, 2 n-12 V 1 12 



(A3) 



167T 



F(ki,k 2 , k 3 , /C4)feifc 2 fc 3 fc 4 



x max[0, min(fci + k2, + fei) 

— max(|fci — /C21, 1^3 — ki\)]dkidk2dk 3 dk4 L . (A5) 



where | A;i — A; 2 1 < Q12 < ki+k 2 - The total number density 
integrated over two momenta can then be found using 



sin 0\2dQ\2 
47rfci27rfc 2 gi2rfgi2 = AitklA-Kkl , (A4) 



This integral is better suited to computational eval- 
uation since it is four-dimensional [rather than the 
six-dimensional Eq. IjAip ], and the term introduced to 
compensate for the angular integral has a relatively 
simple form. 
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